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Abstract 

Using semiclassical approximation method, Bose-Einstein condensation (BEC) of a relativistic 
ideal boson gas (RIBG) with and without antibosons in three-dimensional (3-D) harmonic traps 
is investigated. The BEC transition temperature Tc and the Helmholtz free energy at Tc are 
calculated. The effect of the rest mass of the boson on the properties of the system is also studied. 
We find that Tc of RIBG is higher than that of the nonrelativistic approximation. The RIBG with 
antibosons is also investigated and it is found that the Helmholtz free energy of the system with 
antibosons at Tc is lower than that of the system without antibosons. It implies that the system 
with antibosons is more stable. 
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I. INTRODUCTION 



Bose-Einstein condensation (BEC) has attracted a great deal of attention since its the- 

n 

oretical prediction by Einstein in 1925 Initially, in the theoretical research on the 
phenomenon, the relativistic corrections are neglected near the BEC transition temperature 
Tc 2|-y] . Lately, it is found that in the universe some systems consist of bosons of very small 
rest mass. For example, the rest mass of a pair of neutrinos is about 10 ^"g in the neutrinos 
system For such boson systems, the relativistic effects should be considerable. So 

particular attention has been given to the behavior of the relativistic boson gas. 

The properties of the relativistic ideal boson gas (RIBG) have been discussed in early 
papers |6l-ll0l|. In 1965, Landsberg and Dunning-Davies studied the critical temperature 
and the specific- heat anomaly at Tc of the RIBG [g], Q]. On the basis of the above works, 
Beckmann et al. investigated the behavior of RIBG with and without mass in all space 
dimensions (d) [10]. It was found that for the massless gas, the specific- heat has a gap in 
d > 2, while for the massive gas, there was a gap in c? > 4. Although the research on the BEC 
in the relativistic system has a long history, a complete treatment using the quantum field 
theory has not been carried out until 1980s, in which the pair creation of boson-antiboson 



is considered 
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15|. Haber and Weldon firstly gave the high-temperature expansions for 



thermodynamic functions of the RIBG by taking into account antibosons 



11| . Frota and 



Singh et al. exhibited in detail the thermodynamic properties of d-dimensional system with 



antibosons [12 
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15| . Moreover, the more exact treatments of BEC in RIBG are given 



16l |. Recently, a relativistic BCS-BEC crossover at zero temperature and finite temperature 



17H19|. 



has been studied by the quantum field theory 

It is worth noticing that the previous investigations on the properties of the relativistic 
boson gas were mainly focused on the systems in the absence of external potential. As we 
know, the external potential is very important to control the characteristics of the Bose 



gas [20 



21| . Su et al. discussed the BEC of RIBG in a general external potential by using 



classical Hamiltonian [S] • In this paper, we study the behavior of the RIBG with and without 
antibosons in 3-D harmonic traps using the quantum theory. 

The outline of this paper is as follows. In section [Tll we study the BEC transition 
temperature and Helmholtz free energy at Tc. Meanwhile, we also investigate the effects of 
the different values of mass on the BEC of the RIBG. In Sec. Illlt we introduce antibosons 
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and discuss the BEC in such a system with antibosons, and Sec. IIVI is a brief summary. 



II. RELATIVISTIC BOSE-EINSTEIN CONDENSATION 

Consider a RIBG system of N bosons in 3-D harmonic traps, and the potential function 
V{r) can be expressed as 



(1) 



2 v--^" ■ - y- 

where m and r are the rest mass and coordinate of each boson, respectively. Ux, ooy and 
Uz are the frequencies of anisotropy harmonic traps. By means of the standard Rayleigh- 



Schrodinger perturbation 

by 

in which 
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23|, the total energy of each boson -En^.n^.n^ can be given 



(2) 



^(0) = Smc^ + h 

3h^ " 



(3) 



(4) 



where h is the Planck constant, and c is the speed of light. Ux-, fiy and are 3-D quantum 
numbers of the boson energy. 

To describe RIBG in grand canonical ensemble, the usual expression for the number of 
particles in statistical mechanics is 
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"^a; •X^u ."^z 



(5) 



where P = l/ksT, is the Boltzmann constant, and fi is the chemical potential. p„ is the 
average number of bosons in the state of energy -En^.n^.n^- The number of particles in the 
ground state Nq becomes to take on macroscopic values corresponding to the onset of BEC. 
It is convenient to separate out the lowest energy Eq^q^ from the sum Eq. ([5]), so one can 
write 
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when BEC occurs, the chemical potential equals to the energy of the lowest state -Eo.o.O) i-e., 
fic — ^ 3mc^ + h {cOx + cOy + coz) /2 + 3h^ (col + ujy + cu^) /32mc^. Eq. ([6]) can be reduced as 

where 



= h{nxUx + UyUy + ra^w^), (8) 
3/^2 

^^^^ = Ye™^ ^ ^^^^ ^ "'^'^^ ^ ^""'^^ ^ ^^^^ ^ ■ 

In order to evaluate the above sum explicitly, it is assumed that the energy spacing 
becomes smaller when N ^ oo. So sums ([7]) become integrals as 

N -No= / ^ . (10) 

Jo exp[/3(£(o) + £(i))]-l ^ ^ 

This approximate transformation is good when trapped particles become very large and 
ksT ^ f^x, f^y, fvjJz- So the BEC transition temperature of the RIBG can be obtained by 
taking A^^q — )■ 0, one gets 

f°° duxdriydnz 
N = / — - — fill 

io exp[/3,(e(o)+e«)]-l' ^''^ 

where Pc = ^/ksTc. Fig. 1 displays the exact Tc numerically versus for different values 
of rest mass. It is found that the BEC transition temperature increases with the decrease 
of the rest mass of the particle, which is similar to the case of a uniform Bose gas It 
is explained that the relativistic effect is obvious for the boson system with very small rest 
mass. The transition temperature is compared with the result based on the nonrelativistic 
approximation (see Fig. 2). It is shown that the relativistic effect results in the increasing 
of the BEC transition temperature, which disagrees with the existing result 

The free energy at the critical point can be also investigated. We now discuss the 
Helmholtz free energy at in harmonic traps. According to the partition function of 
the system, the Helmholtz free energy can be expressed as 



poo 

F = ficN + ksTc / dn^dnydn^{\n[l - exp [-/3c(^„„nj„n, - /^c)]]}, (12) 

where /3c is obtained numerically for a given value of and -En^.n^.n^ is given by Eq. (j2]). 
Fig. 3 shows the Helmholtz free energy at Tc as a function of for different values of the 
rest mass. It can be seen that the smaller the rest mass of boson is, the lower the Helmholtz 
free energy is, namely, the system is more stable. 



III. ANTIBOSON AND BOSE-EINSTEIN CONDENSATION 

To our knowledge, at sufficiently high temperature the quantum field theory requires 
consideration of particle-antiparticle pair production. If is the number of antibosons, 
the system is governed by the conservation of the number Q = N — N, rather than of 



the numbers N and separately. According to Ref 221], E = ±En. For bosons we have 



E = +En and for antibosons we have E = —En. So N = J2nPn replaced by 



Q = N-N=Y,iPn-Pn) 



(13) 



E 



Le 



where p„ is the average number of antibosons in the state of energy —En^,ny,nz- Since the 
number of bosons (antibosons) in various states must be positively defined, we have \fi\ < fic- 
Using the similar method as in Sec. HIl the BEC critical temperature can be obtained by 



N-N (14) 

1 1 

exp[/3c(£(°) + e(^y)] - 1 ~ exp[/3c(£(°) + e^-^) + 2fx,)] - l_ ' 

Fig. 4 gives the behavior of Tc numerically obtained from Eq. (fT^ for the different values 
of the rest mass. It can be seen that the behavior of Tc is similar to the case of the system 
without antibosons. However, comparing to the result of the system without antibosons, 
we find that Tc of the system with antibosons is higher (see Fig. 5). This means as the 
decreasing of the temperature, BEC of the system considering antibosons occurs firstly. It 




drixdnydriz 
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implies that the system with antibosons is more stable, i.e., has a lower Helmholtz free 
energy at Tc. This will be shown to be the case indeed. 

The Helmholtz free energy of the system with antibosons can be also given by 

F = fi,{N-N) (15) 

POO 

+ ksTc / dn^dnydn^{\n[l - exp [-(3c{En^^ny,n, - /ic)]] + ln[l - exp [- l3c{En, ^ny,n, + yUc)]]}, 
Jo 

where F is just Eq. f[T^ but with the second log term caused by antibosons. In Fig. 6, we 
find that for a given value of N the smaller the rest mass is, the more stable the system is. 
However, it is found that F of the system considering antibosons is always lower than that 
of the system without antibosons, which proves the above speculation (see Fig. 7). 

IV. SUMMARY 

In this paper, we have studied the BEC of the relativistic Bose gas with and without 
antibosons in harmonic traps by the quantum theory. It can be seen that the BEC transition 
temperature increases with the decrease of the rest mass of the particle. By using the 
quantum energy spectrum, we find the relativistic BEC transition temperature is higher 
than that of the nonrelativistic approximation. Tc of the system with antibosons is higher 
than that of the system without antibosons, and the calculation of the Helmholtz free energy 
at Tc further reveals that the system with antibosons is more stable. 
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Figure Captions 

Fig. 1. The BEC transition temperatures TcS (in units of K) versus the number of 
particles N without antibosons for different mass and the harmonic frequencies ui^ = ~ 
10^^ Hz^ Uz = 2* 10^^ Hz. The sohd, dotted and dashed hnes represent the various values 
of mass m = 3 * 10~^"^ kg, 3 * 10^'^'^ kg and 3 * 10^'^^ kg. respectively. 

Fig. 2. Difference between the BEC transition temperatures T^s (in units of K) for 
relativistic and nonrelativistic approximation with the mass m = 3 * 10~^^ kg, harmonic 
frequencies uix — i^y = 10^^ Hz, uj^ ^ 2* 10^^ Hz. The solid and dotted hnes stand for 
the relationship between the transition temperatures T^s and the number of particles N for 
nonrelativistic and relativistic cases, respectively. 

Fig. 3. The Hclmholtz free energy F (in units of /xiV) versus the number of particles N 
without antibosons at Tc for different mass and the harmonic frequencies ui^ — ijJy — 10^^ 
Hz, Uz — 2 * 10^^ Hz. The sohd, dotted and dashed lines represent the various values of 
mass m = 3 * 10~^^ kg, 1 * 10~^^ kg and 7 * 10""^^ kg, respectively. 

Fig. 4. The BEC transition temperatures TcS (in units of K) versus the number of 
particles N with antibosons for different mass and the harmonic frequencies ujx = <-^y = 10^^ 
Hz,ujz = 2 * 10^^ Hz. The solid, dotted and dashed lines represent the various values of 
mass m = 3 * 10"^*^ kg, 3 * 10~^^ kg and 3 * 10~^^ kg, respectively. 

Fig. 5. Difference between the BEC transition temperatures TcS (in units of K) with 
and without antibosons for the mass m — 3* 10~^^ kg, harmonic frequencies uix — ijJy — 10^^ 
Hz, Uz — 2 * 10^^ Hz. The sohd and dotted hnes stand for the relationship between the 
transition temperatures TcS and the number of particles N for the cases without and with 
antibosons, respectively. 

Fig. 6. The Helmholtz free energy F (in units of iiN) versus the number of particles 
N with antibosons at Tc for different mass and the harmonic frequencies uix — iiJy — 10^^ 
Hz, u!z — 2 * 10^^ Hz. The sohd, dotted and dashed lines represent the various values of 
mass m — 3* 10~^^ kg, 1 * 10~^^ kg and 7 * 10~^^ kg, respectively. 

Fig. 7. Difference between the Helmholtz free energy F (in units of fiN) with and 
without antibosons at T^. for the mass m = 3 * 10^^ kg, harmonic frequencies = oOy = 10^'^ 
Hz, Uz = 2 * 10^^ Hz. The solid and dotted lines stand for the relationship between the 
Helmholtz free energies F and the number of particles for the cases without and with 
antibosons, respectively. 



8 




N 




N 




N 



